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ABSTRACT 



Observed spectra of Active Galactic Nuclei (AGN) and luminous X-ray 
binaries in our Galaxy suggest that both hot (~ lO^K) and cold (~ lO^i^) 
plasma components exist close to the central accreting black hole. Hard 
X-ray component of the spectra is usually explained by Compton upscattering 
of optical/UV photons from optically thick cold plasma by hot electrons. 
Observations also indicate that some of these objects are quite efficient in 
converting gravitational energy of accretion matter into radiation. Existing 
theoretical models have difficulties in explaining the two plasma components 
and high intensity of hard X-rays. Most of the models assume that the hot 
component emerges from the cold one due to some kind of instability, but 
no one offers a satisfactory physical explanation for this. Here we propose a 
solution to these difficulties that reverses what was imagined previously: in 
our model the hot component forms first and afterward it cools down to form 
the cold component. In our model, accretion flow has initially a small angular 
momentum, and thus it has a quasi-spherical geometry at large radii. Close to 
the black hole, the accreting matter is heated up in shocks that form due to 
the action of the centrifugal force. The hot post-shock matter is very efficiently 
cooled down by Comptonization of low energy photons and condensates into 
a thin and cool accretion disk. The thin disk emits the low energy photons 
which cool the hot component. All the properties of our model, in particular 
the existence of hot and cold components, follow from an exact numerical 
solution of standard hydrodynamical equations — we postulate no unknown 
processes operating in the flow. In contrast to the recently discussed ADAF, the 
particular type of accretion flow considered in this Letter is both very hot and 
quite radiatively efficient. 
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1. Introduction 



There is a general agreement that the observed properties of galactic black hole 
candidates and of active galactic nuclei (AGN) could be best explained in the framework 
of accretion disks around black holes. However, no theoretical accretion disk model could 
explain all the basic properties of these sources. In particular, the best known, standard 
Shakura & Sunyaev (1973) disk model (SSD), predicts a temperature of accreted matter far 
too low to explain the hard X-ray emission {hu ^ lOkeV) that is observed. The observed 
hard X-ray emission could be explained by postulating the existence of a very hot plasma, 
with electron temperature Tg ~ 10^ K, in which soft photons emitted by SSD are boosted 
to higher energies by inverse Compton effect. The question is, how does such a hot plasma 
form in black hole accretion flows. 

The very popular disk-corona (DC) model (e.g., Liang & Price 1977; Haardt & 
Maraschi 1993) postulates the existence of the T^. ~ 10*^ K plasma in the form of a hot 
corona above the cold disk. Because the physics of DC models is largely ad hoc, a typical 
specific DC model contains a set of free tunable parameters. 

In the Shapiro, Lightman & Eardley (1976, SLE) hot, optically thin accretion disk 
model, ions are heated by viscous dissipation of their orbital energy, and inefficiently cooled 
by the Coulomb interaction with electrons. Thus, the ions have temperature close to the 
virial temperature, 10^^ — 10^^ K. Electrons are very efficiently cooled by a variety of 
radiative mechanisms, and this reduces their temperature to about 10^ — lO^'' K, which is 
sufficient for explanation of the hard X-ray radiation. The SLE model is, however, violently 
thermally unstable and therefore it cannot describe really existing objects. 

Detailed models of black hole optically thin accretion flows in which cooling is 
dominated by advection (ADAF) have been recently constructed in many papers (see recent 
reviews in Abramowicz, Bjornsson & Pringle 1998 and Kato, Fukue & Mineshige 1998). 
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ADAFs are hot, with the electron temperature about 10^ — 10^° K, and underluminous, 
L <^ LEdd- Here LEdd = 1-3 x 10^^ (M/M0) erg s"-^ is the Eddington luminosity 
corresponding to the black hole mass M . No ADAF solution is possible above a limiting 
accretion rate that is roughly about 0.1 MEdd = 0.1 LEddjf? ■ In some black hole sources, 
however, observations point to accretion rates and radiative efficiencies that are much 
too high for the standard ADAF model to explain (see review in Szuszkiewicz, Malkan & 
Abramowicz 1996). On theoretical side, no satisfactory model of an SSD-ADAF transition 
has been worked out and thus, in order to account for the presence of cold matter together 
with hot ADAF plasma, one still uses phenomenological arguments (see Kato & Nakamura 
1998). 

A model of inhomogeneous inner region of accretion disk was proposed by Krolik (1998). 
In this model the accretion flow consists of clouds moving in a hot, magnetized intercloud 
medium, which can principally explain the emittion of hard X-rays. Such a structure 
could result from a dynamical photon bubble instability in radiation pressure supported 
disks (Spiegel 1977; Gammie 1998). Due to a signiflcant complexity in describing of the 
inhomogeneous medium, the model contains a number of phenomenological assumptions. 
In particularly, the assumptions of stability of such a conflguration is questionable. 

Motivated by the above difficulties of the existent models for black hole accretion flow 
in explaining the co-existence of the hot and cold components, we have constructed a new 
model that has the following properties. (1) Both very hot, hard X-ray emitting gas, and 
sufficiently cool gas, consistent with the detection of the fluorescent iron line, are present 
very close to the central black hole. (2) Signiflcant part (up to ~ 50%) of total luminosity 
(~ 10^"*^ — X^^'^hEdd) of the object is emitted in hard X-rays. (3) Flow is stationary and 
globally stable for large accretion rates (M ~ MEdd)- 

The key element of the model is that accretion matter has initially a low angular 
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momentum. We keep in mind two kind of objects, in which a low angular momentum 
accretion onto black hole may occur. First, black holes, which are fed by accretion from 
a winds blowing of OB star in binary systems (Illarionov & Sunyaev 1975a). The most 
popular object of this kind is X-ray binary Cygnus X-1 (see Liang & Nolan 1984). Second, 
luminous X-ray quasars and AGN, where the central supermassive black hole is fed by the 
matter, which is lost from stars of slowly rotating central stellar cluster (Illarionov 1988). 
In its basic respects, the geometry of our model closely resembles that of the DC model, but 
there are also significant differences. Our models is also very different from that recently 
discussed by Esin (1997) although they both stress the importance of cooling by inverse 
Comptonization. 

2. Low angular momentum accretion 

Let us consider matter with a low characteristic specific angular momentum £, accreted 
quasi-spherically onto the central black hole. By low l we mean that which corresponds 
to the Keplerian orbit with radius tq = II? JTgC^ in range, Sr^ < ''"o ^ lOOr^. Here 
r^ = 3 X 10^ (M/Mq) cm is the gravitational radius of black hole with the mass M. 
Matter with angular momentum smaller than that in the indicated range could not form 
an accretion disk around the Schwarzschild black hole because the innermost stable orbit 
around such a hole is localized at Vms = 3rg. Matter with i higher than in the indicated 
range could form an accretion disk that extends from the black hole to large radial distances, 
r ^ lOOr^: this would correspond to the previously studied SSD. 

At large radii, r ^ tq, the low angular momentum accretion flow closely resembles 
spherical Bondi accretion (Bondi 1952). Approximated models of spherical accretion 
onto a luminous X-ray central source was studied by several authors (e.g. Ostriker et al. 
1976; Bisnovatyi-Kogan & Blinnikov 1980). It was shown by Igumenshchev, Illarionov 
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& Kompaneets (1993) that inside of the Compton radius tq = 10^(10^-ft'/Tc)rg, where 
the Compton temperature Tc is determined by the 'average' photon energy of the source, 
accretion is almost spherical and supersonic. Our model shows that at smaller radii 
r ~ ro, the flow significantly deviates from the spherical accretion. Fluid elements tend to 
cross the equatorial plane at the radius which corresponds to the Keplerian orbit for the 
angular momentum of the element. This leads to formation of shocks above and below the 
equatorial plane at r ^ tq. By crossing the shocks, protons reach the virial temperature 
Tp K, 10^^(rg/r) K. In the presence of soft photons from the thin accretion disk, electrons 
in the post-shock region are efficiently cooled by inverse Comptonization. These cold 
electrons also efficiently cool protons via Coulomb collisions. Such a plasma undergoes a 
runaway cooling due to the intensive bremsstrahlung-Compton processes at layers, where 
the Compton y-parameter reaches the order of unity. At these layers protons loose the most 
of their thermal energy, and the plasma condenses into the thin and cold disk. 

The thin disk spreads to very large radii, r ^ tq, due to the viscous diffusion (von 
Weizsacker 1948, also see Pringle 1981) from the region of condensation at r ~ tq. It 
is convenient to consider two different parts of the thin disk. The inner part, t ^ Tq, 
is an accretion disk, where matter moves inward and angular momentum is transported 
outward. Matter mainly enters the black hole through this accretion disk, which means 
that radiative efficiency of the system is as high as in the standard SSD model (~ 6% for 
the Schwarzschild black hole). At the outer part of the disk, r ^ tq, angular momentum is 
transported outward, removing its excess from the accretion flow to large radii. 

3. Numerical method 

To study details of the model briefly described in the previous Section, we have 
simulated the low angular momentum accretion of weakly magnetized plasma onto black 



hole with the help of two-dimensional time-dependent hydrodynamical calculations. Our 
code is based on the PPM hydrodynamical scheme (Colella & Woodward 1984), and solves 
non-relativistic Navier-Stokes equations in spherical coordinates, assuming the azimuthal 
symmetry. We calculate separately the balance of electrons and protons internal energies 
for the electron-proton plasma. The energies of electrons and protons are coupled through 
the Coulomb collisions. Electrons are cooled by bremsstrahlung and Compton mechanisms. 
The cold and thin accretion disk at the equatorial plane emits soft photons needed for the 
Compton cooling. The thickness of the disk is not resolved in our models. We approximately 
calculate the energy release in the disk, and the correspondent radiation flux, which directly 
affects the efficiency of the Compton cooling, using the standard SSD model. In the model, 
the accretion rate in the cold disk equals to the condensation rate of the hot plasma into 
the disk, which is directly calculated in numerical simulations. 

We neglect the multiple photon scattering processes, when calculating the Compton 
cooling of plasma. This approximation is quite reasonable for the optically thin flows 
(r <^ 1), which are characteristic for our models. We do not solve the transfer equation for 
radiation emitted by the thin disk. However, we calculate the radiation density (which is 
used to find the plasma energy losses) at each point above and below the disk by neglecting 
the absorption of photons. 

In our simulations, protons are heated by shocks, adiabatic compression and viscous 
dissipation. We take into account all the components of the viscous stress tensor 
corresponding to shear in all directions. The bulk viscosity is not considered. The kinematic 
viscosity coefficient is taken in the standard a-parameterization form: u = acl/Qx, where a 
is a constant, c^ is the isothermal sound speed, and D.k = {rg/rY^'^c/\^rg is the Keplerian 
angular velocity. We use a = 0.1 in numerical models. No artificial numerical viscosity was 
used in the calculations. 



-9- 



At the outer boundary we assume a supersonic matter inflow with the spherically 
symmetric distributed density and the specific angular momentum distributed in the 
polar direction consistently with rigid rotation: i{9) = ^max sin^(6'). The parameter imax 
determines (roughly) the maximum radius, r^ = 2if^^^/rgC^, inside which the condensation 
of the hot plasma to the thin disk takes place. For convenience, we will use the parameter r^ 
rather than imax, when describing results of numerical simulations. At the inner boundary, 
Tin = Svg, and at the equatorial plane, where the condensation occurs, we assume total 
absorption of the infiowing matter. 

4. Results and discussions 

We follow the evolution from an initial state, until a time when a stationary fiow 
pattern is established. Models show a strong dependence on two parameters: accretion rate 
M, and r^. We have calculated several models with different M and r^. Two examples 
of stationary accretion flows with r^ = SOr^ and two different M = 0.25 and O.bMEdd 
are presented in Fig. 1. One can clearly see two stationary shock structures, which are 
developed in each model. The structures are more compact for larger M. The nature of 
inner shock, at radial distances r < r^ above and below the equatorial plane, has already 
been discussed in Section 2. In these shocks, the supersonic accretion flow is slowed down 
before condensation to the thin disk. We expect that, in the hot post-shock plasma, the 
inverse Comptonization of soft photons emitted from the thin disk provides the main 
contribution to the hard X-ray luminosity of this type of accretion flow. In this region the 
maximum proton temperature Tp ~ lO^^i^ and the electron temperature Tg ~ lO^i^ in both 
models shown in Fig. 1. The distribution of T^ in the post-shock region is quite uniform. 
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Note, that we do not resolve in our calculations the regions close to the thin disk 
surface, where the main part of energy of the matter being condensed is released, and 
where the Comptonized spectrum is formed. We are limited in space resolution by the 
sizes of numerical cells, but the scale of condensation region is much less than the cell's 
sizes. This lack of resolution can be demonstrated by the estimation of the y-parameter, 
calculated as the integral y = J{kTe/mec'^)dT, where r is the Thomson optical depth and 
the integration is taken in 2;-direction from the equatorial plane till the outer boundary. For 
the larger accretion rate model, M = 0.5MEdd, the maximum value of y is only of order 
of 0.1, whereas the optical depth in vertical direction does not exceed 0.5. The spectrum 
of the escaping radiation is most probably formed at layers with y — 1 and r ~ 1, and its 
hardness is determined by the temperature of these layers. To resolve the layers, where 
the spectrum is formed, one should increase the resolution of numerical scheme at more 
than ten times in comparison with the present one (we currently use rir x ng = 150 x 100). 
As an alternative approach we propose to study the vertical structure of the condensation 
region using a one-dimension numerical scheme with the boundary conditions taken from 
the two-dimension models. 

Radiative shocks of discussed type could be thermally unstable (see Saxton et al. 
1998), and their instabilities may drive quasi-periodic oscillations with time scale of order of 
lOvg/c, which could explain the observed high frequency variability in the X-ray flux (Cui 
et al. 1997). 

The outer shock shapes depend mainly on M and vary from an oblate spheroid to 
a torus (see Fig. 1). The formation of this shock is connected to the presence of the 
centrifugal barrier: supersonically moving matter is stopped at r ;:^ r^ by the centrifugal 
force. The possibility of such kind of centrifugal driving shocks in accretion flows was flrst 
pointed out by Fukue 1987. In the pre-shock region, between the outer boundary and 
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the outer shock, the viscous transport of angular momentum is not important due to the 
supersonic accretion velocities. Slowing down of matter inflow on the outer shock makes 
the viscous transport of angular momentum important in the outer post-shock region. The 
angular momentum is effectively transported outward there. This outward transport, from 
the rapidly rotating inner parts to the slowly rotating outer ones, is balanced by the inward 
advection of angular momentum by the inflowing matter. This balance, as well as the 
efficient cooling of plasma via the Compton mechanism, play a significant role in limiting 
the size of the outer post-shock region and creation of the stationary shock structures. We 
found, that in the case when the Compton cooling is (artificially) switched off, no stationary 
solutions of this type are possible. 

The dependence of the flow structure on the parameters r^ and M can be used for 
explanation of the observed spectral variability of galactic black hole candidates. For 
example, the change of spectral state from 'hard' to 'soft' in the case of Cygnus X-1 (see 
Tanaka & Lewin 1995) can be explained in terms of variation of ra- In the hard state (small 
Td) most of the accreting matter goes through the inner shocks at radius ~ lOr^. In this 
case a significant part of radiation comes in hard X-rays. In the soft state (large r^^) the 
accretion matter comes to the thin disk at radii ^ lOOr^. Close to the black hole, accretion 
through a thin disk takes place only. Obviously, hard X-ray emittion is suppressed in this 
case. Both parameters r^ and M are very sensitive to conditions of mass exchange between 
binary companions in the case of a wind fed accretion (Illarionov & Sunyaev 1975b). A 
small change in the velocity of the wind of an OB star companion leads to a significant 
variation of both r^ and M with amplitudes that could explain observed variability. 
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5. Conclusion 



New model of accreting black holes have been proposed, which self-consistently explains 
the existence of two-component (hot and cold) plasma in the vicinity of black holes. The 
model assumes that at large radial distances the accretion flow has a quasi-spherical 
geometry and a low characteristic angular momentum. The main parameters of the model 
are the accretion rate and the amount of angular momentum carrying by accretion flow. 
The model can naturally explain, without phenomenological assumptions, the origin of hard 
X-ray excess observed from the stellar mass black holes (Galactic black hole candidates), 
as well as the supermassive black holes (quasars and AGN). More hydrodynamical and 
radiative transfer simulations are needed to construct detailed spectra of accreting black 
holes, which will be compared with observations. 

This study was supported in part by the grant from the Royal Swedish Academy of 
Sciences, RFBR grant 97-02-16975 and NORDITA under the Nordic Project 'Non-linear 
phenomena in accretion disks around black holes'. We would like to thank the referee for 
many helpful suggestions. 
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Fig. 1. — Low angular momentum accretion onto black hole. Distributions of density p 
and velocity vector field are shown in the meridional cross-section for two models with 
M = 0.25MEdd (left) and 0.5MEdd (right). Note the different scaling in the axes of the left 
and right panels. Vertical axis coincides with the axis of rotation. Black hole locates in the 
origin (0,0). The contour lines are spaced with Alogp = 0.1. Cold thin disk, which is not 
resolved in our numerical models, locates at the equatorial plane, z = 0. Condensation of 
hot matter to the cold disk takes place inside the radius ^ r^ = SOr^ at the equatorial plane. 
Accretion flow forms a two-shock structure (shocks are seen as concentrations of density 
contours). The structure is more compact for the higher accretion rate model (right panel). 
In the inner shocks, the supersonic accretion flows, which locate close to the polar directions, 
are slowed down before condensation to the thin disk. The outer shock originates due to 
a centrifugal barrier, where the infalling matter is stopped at the equatorial region by the 
centrifugal force. 



